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A dc current can be pumped through an interacting system by periodically varying two inde- 
pendent parameters such as magnetic field and a gate potential. We present a formula for the 
adiabatic pumping current in general interacting systems, in terms of instantaneous properties of 
the system, and find the limits for its applicability. This formula generalizes the scattering approach 
for noninteracting pumps. We study the pumped spin in a system that exhibits the two-channel 
Kondo effect as an application of the adiabatic pumping formula. We find that a quantized spin 
of h is transferred between the two channels as the temperature approaches zero, and discuss the 
non-Fermi liquid features of this system at finite temperatures. 

PACS numbers: 72.25.-b, 73.23.-b 



Introduction and Conclusions. — The scattering ap- 
proach of Brouwer Q for pumping through a finite, pos- 
sibly disordered region of noninteracting electrons, which 
followed a work by Biittiker, Pretre and Thomas 
(BPT), tremendously enhanced the understanding of 
time dependent transport in mesoscopic systems. The 
derivation of the Brouwer formula Q and its applications 
are well established 0, LU , and recent studies considered 
special cases of interactions in quantum dots H@ and in 
Luttinger liquids However, an apt formulation of 

the pumped current in the general case when interactions 
between the electrons are involved (beyond the Hartree 
level) has been lacking. 

In this Letter we develop a generic formula [Eq. J2J] 
which expresses the pumped current through a region of 
interacting electrons in the adiabatic limit. The pumped 
current is expressed in terms of an instantaneous lin- 
ear response function, calculated at every step along the 
pumping trajectory. When the motion along the tra- 
jectory is sufficiently slow, the current can be found by 
integrating the contributions of each step [see Eq.Q]. 
We also rewrite Eq. (|2J for the case of a quantum dot 
connected to noninteracting leads in terms of dot prop- 
erties [Eq. from which the noninteracting S-matrix 
formula |2j follows as a special case. 

Finally, we consider as an example spin pumping in 
the vicinity of the two-channel Kondo (2CK) fixed point 
which is perturbed by two pumping parameters, mag- 
netic field and channel anisotropy. Due to the non Fermi 
liquid (NFL) nature of the 2CK fixed point, the stan- 
dard scattering approach can not be applied straightfor- 
wardly. Using Eq. (0) we calculate the spin pumped from 
one channel into the other, as function of the parametric 
trajectory and temperature T. We find that as T — ► 
the pumped spin is quantized in units of h per period for 
trajectories which surround the NFL fixed point [para- 
graph following Eq. (0] , while at finite temperatures this 
quantization is not accurate [Eq. Jj5J], At finite tempera- 
ture, when the trajectory is sufficiently close to the fixed 
point, the temperature dependence of the pumped spin 
reflects the NFL physics of the 2CK [Eq. JTUJ)]. 



Generic adiabatic pumping formula. — We analyze a 
mesoscopic conductor which is described by an Hamil- 
tonian Hx 1 ,x 2 ,... that depends at least on two external 
parameters X\ (t) and X 2 {t) which are varied periodically 
and slowly in time. These parameters can be for exam- 
ple two metallic gates and/or an external magnetic field. 
By assumption the conductor may be divided into left 
(L) and right (R) contacts whose Hamiltonian does not 
depend on the parameters, and a central region whose 
Hamiltonian does depend on the parameters. Interac- 
tions may take place everywhere in the conductor. 

In a pumping cycle the set of the parameters X(t) = 
(X 1 (t), X 2 {t), . . . ) is varied periodically in time and de- 
fines a closed trajectory, C, in the parameter space. The 
change in the parameters may produce a current Jj in 
contact- 7 (j = L,R) |9j. To calculate Jj at a point Xq 
along the trajectory due to a small and slow change, 8X, 
of one of the parameters, we assume that SX = SX^e . 

The current Jj(fl) is then found by a linear response 
perturbation theory with respect to the infinitesimal per- 
turbation H' = jJ^8Xne inT . Including the possibility 
that the current operator, Jj, may depend explicitly on 
X, we obtain the total response to linear order in SXn'- 
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with = fi — i0 + . Since Eq. Q for the current is 
a first order expansion in the harmonic perturbation, 
the time evolution should be understood as 0(t) = 
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formed with (f2 independent) instantaneous eigenstates 

The charge SQ(j, fi) = Jj(£l)/ (itt) entering the central 
region through contact j is SQ(j, f2) = e dr ^ 
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, the emissivity into contact j, is given by: 
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The last equality follows since in the static limit, f2 = 0, 
no current flows through the contacts. This yields 0, 
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Notice that while here we take the limit O — > 0, in prac- 
tice the emissivity does not depend on Q for f2 < l/r r , 
where r r is a characteristic relaxation time. The depen- 
dence of the emissivity on time is only through the loca- 
tion of Xq on the trajectory C 

Using Eq. f° r the emissivity related to each param- 
eter, the charge pumped per period, corresponding to a 
trajectory C in the parameter space, is given by |l| 
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where S is the area bounded by C and the effective "vec- 
tor potential" and "magnetic field" are A(j) = e and 

B(j) = (V x A(j)) 3 = - respectively. 

To determine the validity regime of Eqs. J5J) and (JSJ 
we divide the trajectory C into elements of length SX. 
The length, SX, should be smaller than both the radius 

U 3 , ~ 

of curvature along the trajectory, r c = X I X x X 

and r 2 = jx^Xi j d x k dx ^kXj , the length at which 
the second order term in SX is comparable to the first or- 
der term. Here X% — X{ / X , and the Einstein summa- 
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tion convention is understood. In addition, SX = 

and St must be longer than r r , so that after each step 
the system relaxes to a new equilibrium position de- 
termined by the Hamiltonian with the new parameters 
X + dX/dt St. Combining both requirements we find 
that for every X 6 C, the length SX has to satisfy: 
St = SX <C min{r2,r c }. Thus, formu- 
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For a circular trajectory: X(t) = 7*0 (cosfioi, sinOo*)) 



one easily finds that r c = ro, dX jdt = f^o an d the 

condition becomes: f^o <^ t - min{r '2/ro, !}■ 

Adiabatic pumping formula for quantum dots. — 
Next we consider a quantum dot coupled to non- 
interacting contacts, described by a parametric de- 
pendent version of the Hamiltonian discussed in 
Ref. [13, = E k , aeL , R e ka c{ a c ka + Hf({dld n }) + 

J2k, a <£L,R Vka,n(X)cl a d n + h.c. Here c\ a creates an elec- 
tron with momentum fc in channel a belonging to contact 
j, and {d^} form a complete, orthonormal set of single- 
electron creation operators in the dot. In this case the 
current operator is Jj = ^ J2k,aej v ka,n{X)c\ a d n + h.c. 
A straight forward calculation shows that using Eq. (J2J 
the emissivity can be written as {h = 1): 
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where u = h- t[, T = t -^. 
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Here G^ v b (£1, i^, O) is the Keld ysh Green function related 
to the operators a,b — d n , Cka [13J. We note that the de- 
pendance of the Green function on SXq is through H'[t) 
defined above Eq. QJ. 

Since the contacts are noninteracting, the summation 
over k in Eq. |J2J can be carried out [12J, and the emis- 
sivity can be written in terms of SxG r ^ 1 m only. Defining 
a vertex function 

S x G(e, il) = G(e - ~)A x (e, Q)G(e + -), 

where we used matrix notation both for the dot indices 
(bold letters) and for the Keldysh indices (hat, *), with 
A / G r G< 
G a 
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The arguments of /(e), G(e), T(e), and Ax(e,fl) were 
suppressed. nj = J2 a ej,k I du>f(u) [ - ^lmG ka ,ka {u>)] 
is the equilibrium occupancy of contact j calcu- 
lated in the presence of the dot, r£_ m (e) = 
27r Sa6j K,a(e)Pa(e)V r a , m (e), p a (e) is the bare density of 



V ka ,n for e = The 



P Bm ^ is antihermitian. 



states in channel a and V atn {e) 

matrix B x T^ m = ^T 3 n , 

In the noninteracting case one can show that expres- 
sion JSJ) without the first term, — ^y, reduces to BPTs' 
result. The BPTs' emissivity contains the explicit deriva- 
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tive contribution as they calculate the current 

deep inside the reservoir, while in this Letter the cur- 
rent defined above Eq. |j2J is calculated at the entrance 
to the dot. An explicit derivative does not influence the 
pumping charge per period, and formally corresponds to 
a gauge transformation in the vector potential defined 
after Eq. ©. Notice that for an infinite flat band and 
energy independent tunneling couplings = even 
when interactions in the dot are included [l5| . 

Pumping in the two channel Kondo effect. — To demon- 
strate the new features of Eq. @ which includes inter- 
actions, we study a specific example of pumping in a 
2CK system at the exactly solvable Emery-Kivclson (EK) 
line The peculiarity of a symmetric 2CK problem 

is in its NFL behavior at low temperatures 0|- I* 1 the 
presence of external magnetic field, B, the Hamiltonian 
of the 2CK model is 

H 2CK = ^2 e k c\ aj c kaj + J2 JjxSx ■ ^ + gti B HS z . 

kaj j,\ 

The index j = 1,2 represents two channels, S is the im- 
purity spin-l/2-operator, is the spin density in channel 
j in direction A = x,y, z near the localized spin. 

In the present context we consider spin pumping from 
channel 1 into channel 2 by calculating the spin-flavor 
(sf) emissivity = f - using Eq. @ 

with Jj replaced by J s f = ~e^f- = ^ [n 8 f, H 2GK ], where 
n s { = s z (l) — Sz(2) is the spin difference between the 
channels. At the EK line J7i 2 = = 2ttHvf the spin 
charge and flavor sectors commute with the spin-flavor 
(sf) sector that determines the evolution of J s f . 

In the presence of channel anisotropy in the spin flip 
processes, J\y ^ the Hamiltonian of the SF sector 

in a Nambu notation: ^} = tpi) , i = d, sf takes the 

form of a Majorana resonance level (MRL) model 

+\htf d T z ^ d + *t f (0)F% d + ^t>* sf (0)), 

where V = yjT/{2np) (cos(6/2)t z + isin(6/2)Ty), V = 
Ti + T2, r 1( - 2 ) = pJi(2)±_ l^ a i a can °e considered as 
the lattice spacing, h = 2gp B H/T, p = (2nhv F )~ 1 
is the density of states of the chiral fermion field i/> s f, 
and r are the pauli matrices. The operator that de- 
scribes transition of spin between the channels is n s f = 
dxipl { (x)i(j s f(x) |l8| and tpd is a local fermion opera- 
tor. 

The unique properties of the 2CK system can be seen 
by taking 9 = tt/2 and h = in the MRL model. In 
order to describe pumping in the vicinity of this point 
we choose as pumping parameters X\ = h and X2 = 
A = (Ti -r 2 )/r = cos0. 



Using Eq. (J2J) one finds that the spin flavor emissivity 
is given by A = = - de^a(e), where 




FIG. 1: The effective magnetic field B in the (h, A) param- 
eter space. For T « T the effective magnetic field has a 
peak of weight h at (h, A) = (0, 0) and the shape of the peak 
is temperature independent if plotted against the scaled pa- 
rameters O/v'SttT/T, A/y/8-TrT/r). Each of the pumping 
trajectories Ci : (h/ho) 2 + (A/A0) 2 = 1 (full line) and £2 : 
{(h + 3h )/h ) 2 + (A/A0) 2 = 1 (dashed line) is plotted for 
Ti <C TAq and for T-i 3> TA 2 , in the scaled parameter space, 

T ■ 

leading to four different curves C i 3 , i,j = 1,2 (ail trajectories 
are anticfockwise). Denoting the area bounded by C i 3 as S i 3 , 
we see that while S^ 1 contains the entire area of the peak of 
B, correspoding to pumping a spin of exactly h from chan- 
nel 1 to 2, the intersection of the area of the peak with Sj 1 
is empty. On the other hand, the relative area of the peak 

contained in is approximately A °jyp /,r , corresponding 

to the pumped spin ~ SAo / ^JT/Y. The anomalous power is 
a manifestation of the NFL behavior. 

At T = Eq. © gives A = (which is ^ in 

polar coordinates) - the vector potential corresponding 
to an effective magnetic field B = H S 2 (h, A) perpendic- 
ular to the plane of the trajectory. Thus, in a pumping 
period encircling the NFL point (h, A) = (0, 0) anticlock- 
wise, a total spin of exactly h is transferred from channel 
1 to channel 2 0] . The magnetic field B can be obtained 
analytically for any finite temperature. We discuss below 
the different regions of B as a function of T. 

For T/T <C 1, namely for temperatures smaller than 
the Kondo scale, the pumping vector potential can be 
approximated by a(e) = which after Integra- 

tion over e gives A = ifti (g + g^fOi where ipi is the 
trigamma function. (We have neglected the third term in 
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Eq. JJjJ which is curl-less.) The magnetic field B = V x A 
is 

B(r) = ^r/(8^r)F! (r/v^T/r) , (7) 

with the function Fi(x) = ^-^{y^ii^ + y))\ y =x 2 satis- 
fying Fi(0) = 7r/2 and Fi(x) x>1 > q^e- This means 
that the effective magnetic field of strength ~ hT/(16T) 
is concentrated in a circle of radius ~ y/T/T, and decays 

strongly as ^ T ffl , as depicted in Fig. ^ 

At T ~ T the peak size approaches unity and ceases to 
be circularly symmetric due to the second term in Eq. JfjJ). 

For T ^ F the effective magnetic field becomes prac- 
tically independent of — 1 < A < 1, and given by 

%t f t\ 
B ( h > T ) = wT F2 { h Tj> (8) 

with F 2 (x) JLRe^i-^)]. Since 

J_ dxF%(x) = 8, the weight of the pumping peak 
is again H, however the peak is very wide ~ T/T 3> 1. 

The total pumped spin is obtained by performing the 
integral sf_^ 2 = Is d 2 rB(r), where S is the area con- 
tained in the trajectory C. We can easily estimate the 
temperature dependance of the pumped spin, using the 
structure of B(r) described in Eqs. Q and JSJ). 

Consider for example, C\, an elliptic pumping trajec- 
tory (h/h ) 2 + (A/A ) 2 = 1 where A <C 1 < h (see 
Fig. nj . At T — it encircles the origin and therefore 
sfi +2 (T = 0) = h. At low temperatures y/T/T < A , 
taking into account the tail of the peak of B, we obtain 

This is expected from the FL behavior along the trajec- 
tory (no anomalous exponents appear). At higher tem- 
peratures when Ao <C y/T/T <C 1 wc find 

sf^(T) = c hy/T/fA , (10) 

with c of order unity. Here the anomalous exponents 
(~ T~ x / 2 ) of the NFL point become apparent. In Fig.Q] 
the magnetic field is plotted in terms of scaled param- 
eters (h/y/8irT/T, A/y/8nT/T). Using this scaled pa- 
rameters, as long as T < T, the shape of the peak is 
temperature independent, however each trajectory £ ac- 
quires a temperature dependence C T . Let u s consider 
the temperatures T\ and T 2 satisfying y/T\/T <C Ao and 
Ao <C y/T-z/T <c 1 respectively We see that the area 
bounded by C^ 1 contains the entire peak of B while the 
area bounded by Cj 2 contains only a one dimensional 
cut of the peak, (whose radius scales as T/T for the bare 
parameters) explaining the anomalous behavior for T 2 . 

For 1 < T/T < h wc find sf^ 2 (T) = hA . Finally 
for very large temperatures 1 -C ho <C T/T we have 



s f^ 2 ( T ) = ft7r/ioA r/(16T). Such area law [ 0(irh Ao)] 
is expected in this regime, since B is practically constant 
for h ,A < T/T. 
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